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Abstract 



The A'^-dimensional Hamiltonian 



is shown to be quasi-maximally superintegrable for any choice of the functions / 
and lA. This result is proven by making use of the underlying s[(2, R)-coalgebra 
symmetry of Ti in order to obtain a set of (2A^— 3) functionally independent integrals 
of the motion, that are explicitly given. Such constants of the motion are "universal" 
since all of them are independent of both / and lA. This Hamiltonian describes the 
motion of a particle on any A^D spherically symmetric curved space (whose metric 
is specified by /) under the action of an arbitrary central potental and includes 
simultaneously a monopole-type contribution together with centrifugal terms that 
break the spherical symmetry. Moreover, we show that two appropriate choices for 
lA provide the "intrinsic" oscillator and the KC potentials on these curved manifolds. 
As a byproduct, the MIC-Kepler, the Taub-NUT and the so called multifold Kepler 
systems are shown to belong to this class of superintegrable Hamiltonians, and 
new generalizations thereof are obtained. The KC and oscillator potentials on A'^- 
dimensional generalizations of the four Darboux surfaces are discussed as well. 
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1 Introduction 



An A^-dimensional (A^D) Hamiltonian H is said to be (Liouville) integrable [1] if it admits 
the maximum number {N — 1) of functionally independent and Poisson-commuting global 
first integrals. Similarly, an integrable Hamiltonian Ti. is called superintegrable if there 
exists an additional set of k independent constants of motion. It is well-known that 
k < N — 1, and in the case where this bound is saturated {k = N — 1), the Hamiltonian 
Ti is called maximally superintegrable (hereafter MS). When k — N — 2, we shall say that 
Ti. is quasi-maximally superintegrable (QMS) [2, 3] and the Hamiltonian H. is only "one 
integral away" from being MS. 

Superintegrable systems have been thoroughly studied because of their significant con- 
nections with generalized symmetries [4, 5], isochronous potentials [6, 7], and separability 
of the associated Hamilton- Jacobi and Schrodinger equations [8, 9, 10]. When one con- 
siders a Hamiltonian system in arbitrary dimension N and all the integrals of motion 
are imposed to be quadratic in the momenta, the list of known NY) MS systems becomes 
strikingly short. To the best of our knowledge, the only known instances are: 

• The free motion on the simply connected spaces of constant curvature (see, e.g., [11, 
12] for a unified approach in terms of the curvature) which is given by the geodesic 
flow on these spaces. 

• The generalized Kepler-Coulomb (KC) system on the simply connected spaces of 
constant curvature [12, 13, 14, 15, 16], which is a superposition of the A^D KC 
potential with N — 1 "centrifugal" terms. 

• The Smorodinsky-Winternitz system on the simply connected spaces of constant 
curvature [11, 12, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25], that is, the NB harmonic 
oscillator potential on such spaces together with N centrifugal terms. 

• The free motion and the Smorodinsky-Winternitz system on the Darboux space 
III [26], which is the superposition of an "intrinsic" oscillator and A^ centrifugal 
potentials on an A^D generalization of the Darboux surface of type HI, which is a 
Riemannian manifold of nonconstant curvature. 

Note that the above MS systems on spaces of constant curvature do not only include the 
fiat A^D Euclidean case, but also the spherical, hyperbolic, Minkowskian and (anti-)de 
Sitter spaces (for the relativistic ones see [12, 16]). 

On the other hand, if the condition that all the integrals of the motion have to be 
quadratic in momenta is suppressed, the list of A"D MS systems is enlarged with the 
rational and hyperbohc Calogero-Sutherland-Moser models of type Ajv-i [27, 28, 29, 30], 
the nonisotropic oscillator with rational frequencies [31], the nonperiodic Toda lattice [32, 
33], the Benenti systems [34, 35] and, very recently, the KC potential on the 3D Euclidean 
space plus three centrifugal terms [36]. 
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From this viewpoint it is quite natural to wonder whether the number of A^D QMS 
systems is much larger than the number of MS ones (obviously, any MS system is by 
construction a QMS one). The answer to this question has been recently answered af- 
firmatively in [37] as follows: geodesic flows on any ND spherically symmetric (curved) 
space define always a QMS Hamiltonian. Moreover, all these systems present the same 
"universal" and explicit set of {2N — 3) integrals in involution with the Hamiltonian that 
are quadratic in the momenta. The proof of this result is based on the common 5[(2, R) 
Poisson coalgebra symmetry [38, 39, 40] of all these systems, a fact that was also shown 
for the quadratically MS systems on Riemannian spaces hsted above [2, 3, 26, 37, 38]. 

The aim of this paper is to use the same sl{2, M)-coalgebra symmetry in order to 
show that such infinite family of QMS free systems in arbitrary dimension can be further 
enlarged by including the Hamiltonians 



where / and U are smooth functions and j2,bi, . . . ,b]\f are real constants. The proof that 
this Hamiltonian is indeed QMS is presented in section 2. We stress that, despite that 
the spherical symmetry is broken by the "centrifugal terms" bi/qf, the overall 5l(2,M)- 
coalgebra symmetry holds and, as a consequence, the same integrability properties of the 
geodesic flows arc preserved when the potentials are considered. Next sections will be 
devoted to the study of certain particular cases of (1.1), both from the physical and from 
the mathematical point of view. 

Firstly, note that the term |p^//(|q|)^ is just the kinetic energy that defines the 
geodesic flow on a spherically symmetric space, which is generically of nonconstant cur- 
vature (/ playing the role of a "conformal factor" of the flat metric). These are just the 
aforementioned geodesic flows on NY) spherically symmetric spaces [37]. In particular, 
suitable choices for the function / allow us to recover the three classical Riemannian spaces 
of constant curvature [41], A^D generalizations [26, 37] of the Darboux surfaces [42, 43] 
and the Iwai-Katayama spaces [44] (which are themselves a generalization of the Taub- 
NUT metric [45]). All these QMS geodesic flows are fully described in section 3 by using 
the s[(2, R)-coalgebra framework, where as a new result the complete integrability and 
separability of the free Hamiltonian (for any /) is explicitly shown by using spherical 
coordinates. 

Secondly, section 4 is devoted to the study of the central potential U and to the 
analysis of its separability properties which hinges on its Lie-Poisson coalgebra sym- 
metry. In section 5, the A^D versions of the potential U defining the "intrinsic" KC 
and oscillator potentials on generic spherically symmetric spaces is presented by extend- 
ing to arbitrary dimension (sec [26]) the definitions of such potentials on 3D Rieman- 
nian spaces previously introduced in [46, 47, 48, 49, 50]. Moreover, by taking into ac- 
count that the ^^/(2/(|q|)^q^) contribution can be interpreted as a Dirac monopole-type 
term [44, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61] and by adding A^ centrifugal terms to 




(1.1) 
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the KC/oscillator potential new generalized KC and oscillator systems on arbitrary 
spherically symmetric spaces are introduced. These general results are exphcitly illus- 
trated in section 6 by studying in detail the resulting potentials for some specific spaces. 
In particular, known superintegrability results for the MIC-Kcplcr [44, 56, 60, 62, 63, 64] 
and the Taub-NUT systems [44, 59, 65, 66, 67, 68, 69, 70, 71, 72] are recovered within our 
framework, and we also derive further results for new Hamiltonians such as the proper 
KC and oscillator potentials on the A^D Darboux spaces introduced in [37]. This way new 
QMS systems on A^D curved spaces are presented and analysed under a unified integra- 
bility setting, showing that the integrability properties of several well-known systems rely 
on a common sl{2, M)-coalgebra symmetry. 



2 The Hamiltonian H is QMS 



The following result holds: 

Theorem 1. Let {q, p} = {{qi, ■ ■ ■ , Qn), (pi, ■ ■ ■ ,Pn)} be N pairs of canonical variables. 
The ND Hamiltonian 



where f andlA are smooth functions and /i, 6i, . . . , fejv 'J'^e real constants, is quasi-maximally 
superintegrable. The {2N — 3) functionally independent and "universal" integrals of the 
motion for 7i are given by 

i2 ^2^ 



i<i<j ^ \ y J i=i 

N-m+l<i<j ^ \ / J j=Ar_m+l 

where m = 2, . . . , iV and C^^^ — C(^n)- Moreover, the sets of N functions {H, C^"^^} and 
{7i, C(^rn)} {m — 2, . . . , N) are in involution. 

Proof. Let us define the functions 

J_ = q^ J3 = q-P, J+ = p^ + 5]4, (2.3) 

i=i ^» 

where 

N N N 

i=l i=l i=l 
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It is immediate to check that these three functions close the 5[(2,M) Poisson coalgebra 



{J3,J+} = 2J+, { J3, J_} = -2J_, {J_,J+} = 4J3. (2.4) 
Moreover, V, can be written as the following function of the 5l(2, R)-coalgebra generators: 

As a consequence of this coalgebra symmetry, Ti. will Poisson-commute with all the in- 
tegrals C*^"*^ and C{rn) that respectively come from the left and right mth-coproducts of 
the Casimir function for sl(2,]R). This coalgebra symmetry also ensures that the sets of 
TV functions {7i,C(™)} and {7i,C(m)} (1 < m < TV) are in involution (see [3, 38, 40] for 
explicit proofs of all these statements) . 

In this respect, we recall that any Hamiltonian 

n^n{j-,j+,J3), (2.6) 

defined as a smooth function on the sl(2,M) generators is always QMS, and has the 
same set of "universal" integrals (2.2) that are obtained through the mD symplectic 
realization of the Casimir function of sl(2,M). This underlying coalgebraic structure can 
be interpreted as a generalization of the spherical symmetry, which is recovered when all 

Finally, note that the centrifugal terms bi/qf come from the fact that sl{2, R) allows us 
to add such a contribution in the corresponding symplectic realization of the J+ generator. 
Therefore, from an algebraic viewpoint these centrifugal contributions are directly related 
to the kinetic energy term. 



3 Superintegrability of geodesic flows 

The metric of any A^D spherically symmetric space M. can be written as 

ds' = /(|q|)Mq2 = firfidr' + rW), (3.7) 
which is expressed in terms of the following coordinate systems: 

• Generic coordinates q = (gi, . . . , q^) (with dq^ = ^^j^ dg^^), which in our approach 
we shall identify with those appearing in the symplectic representation (2.3). 

• Spherical coordinates with a radial- type variable r = |q| G M'^ (which does not 
usually coincide with the geodesic distance) and — 1 ordinary angular variables 
6j e [0, 2tt) {j = 1, . . . , N — 1) with dQ^ being the standard metric on the unit 
(A^ - 1)D sphere 

df)' = ^d^ljjsin^^fc. (3.8) 

j=i k=i 
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In both coordinate systems /(|q|) = f{r) is any smooth function, usually interpreted 
as a conformal factor of the Euchdean metric ds^ = dq^. We stress that only in the 
Euclidean case q can be interpreted as Cartesian coordinates. The relations between 
both sets of coordinates read 

i-i N-1 
Qj — r cos 9j Y\ sin 9^, 1 < j < N, = r JJ^ sin 9k, (3.9) 

fe=l k=l 

where hereafter any product YlT ^^ch that Z > m is assumed to be equal to 1. The scalar 
curvature of (3.7) turns out to be 

, _ 1) + 2(iV - l)r-7-(r) + (iV - 2)/-(r)^ _ ^ 

where /'(r) = d//dr and /"(r) = d^//dr^. Therefore, we are indeed dealing with ND 
spaces with nonconstant curvature. 



3.1 Free Hamiltonians 

The metric (3.7) provides the free Lagrangian, which characterizes the geodesic motion 
of a particle on in terms of the velocities q or (r, 9j): 

r = ^/(|q|) V = ^firf + E ^1 11 " (^-H) 

\ j=i k=i J 

Then the canonical momenta p,Pr,P9j, conjugate to q,r,9j, can be obtained through a 
Legendre transformation (p^ = dT/dqi) 

i-i 

P = /(|q|)'q, Pr^firff, pe, ^ f{r)'r%l[sm' 9^. (3.12) 

k=l 

Prom (3.9) and (3.12) the relations between both sets of momenta p and Pr^POj are found 
to be 

^ ■ n n cos9.j ^ nCLi sin 9k ^ sin 9^ 

sm 9k cos 9j Pr + }^ cos 9i pe, - i Pe, , 

k=l ^ 1=1 llm=lSinC^m ^^llfc=lSm^^fc 

N-1 ^ N-1 T-r^^-l • O 

sin9kPr + - > , r=Tzi — ——cos9ipei, (3.13) 



k=l 



^ 1=1 nm=lSin^r, 



where 1 < j < N and from now on any sum such that Z > m is assumed to be zero. 
Both sets of phase spaces (q, p) and (r, 9j;pr,P0j) are canonical coordinates and momenta 
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with respect to the usual Poisson bracket 



_dF^dG__dG_dF^^y. fdF^dG__dG_dF\ 
dr dpr dr dpr \ 06 j dpe^ dOj dpg^ J 

This imphes that the geodesic motion is described through the kinetic energy Hamil- 
tonian 

p2 p2 _|_ j^-2]^2 

where is the total "angular momentum" , namely, 

N-l j-l ^ 

>^^=E^in-^- (3-16) 

Therefore the kinetic term on M. arises as an 5((2, ]R)-coalgebra Hamiltonian (2.6): 

provided that all bi = 0. Therefore, this geodesic flow is always QMS for any function / 
and its constants of motion are just (2.2). 

We remark that in the case of free motion we recover the well-known so{N) Lie-Poisson 
symmetry of spherically symmetric spaces. Explicitly, the functions Jy = QiPj — QjPi with 
i < j and i,j — l,...,N span an so{N) Lie-Poisson algebra 

{Jijj Jik} Jjky {'^ijy "^jk} Jikj {^iky Jjk} "^ijy i <i j ^ k. (3.18) 

Hence the integrals (2.2) can be expressed as sums of such "angular momentum" compo- 
nents 

m N 

CM = = E 4- (3-19) 

l<i<i N-m+l<i<j 

Notice that C*-^-* = C(Ar) is the quadratic Casimir of 5o{N) and, in fact, each CM (^or 
C(m)) is the quadratic Casimir of certain rotation subalgebra 50 (m) C so{N). 

By taking into account (3.9) and (3.13) we obtain the generators written in terms 



7 



of the spherical phase space variables; namely, 

■ n /I TT • /I COS sin 

Jij = sm t'i cos t/j sin6^fcpe. - . , Pe. 

i-k n{-Li sin Ok 

+ cos 6li cos > ;_i cosOipoi, 1 <i < j < N, 

Ff™ sin 6*^, 

Ji^ = sin TT sin p0, + cos 6^ V ^^^^^^ ^ cos Oi pe, . (3.20) 

I" 7"^ L, sin 9jn 

k=i 1=1 A ■l-m=« 

By using these expression one can readily write down the constants of motion (3.19) in 
this set of canonical variables. The resulting formulas for C{rn) adopt the very compact 
form 

iV-l j-l 



E 4 n -ir. (3-21) 



j=N-m+l k=N-m+l 

SO that C(jv) = is the total "angular momentum" (3.16). Furthermore, the complete 
integrability determined by the set of N functions {T, C(m)} {fn = 2, . . . , A^) leads to 
a separable system as a set of N equations each of them depending only on a pair of 
canonical quantities. In particular we obtain a common set of — 1 angular equations 
for the free Hamiltonian (3.15) 

C(2)(6'iv-i,p0jv-i) =pL-i' 

C 

Cq){0n-i+i,P9^_,+,) = Pe;v-z+i + -2!!"'^ 1 l^3,...,N-l, 

C(iv)(^l,P,J=pU^^^L^ (3.22) 

sm Oi 

together with a single radial equation corresponding to the Hamiltonian itself T{r,Pr)- 



3.2 Some examples of QMS geodesic flows 

The s[(2, ]R)-coalgebra setting describes all spherically symmetric spaces and, consequently, 
many geodesic motion Hamiltonians that have been studied in the literature by following 
different procedures can be embedded into this integrability framework. In what follows 
we shall present three families of examples whose explicit metric and free Hamiltonian 
are displayed in table 1. 

• The three simply connected Riemannian spaces of constant curvature. If we take 
/(|q|) = 2/(1 + nc^) we recover the spherical (k > 0), Euclidean (k = 0) and 
hyperbolic {k < 0) spaces of constant sectional curvature k (the scalar one is N{N — 
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Table 1: Geodesic flow Hamfltonians with sl{2, R)-coalgcbra symmetry for some particular 
A'"D spherically symmetric spaces in the generic coordinates (q, p) . 



Space Metric ds^ Hamiltonian T = T(J_, J+) 



Euclidean 


dq2 


1 T 


Spherical 


4dq2 
(l+q2)2 

4dq2 
(1 - q2)2 


i(i + J_)V+ 


Hyperbolic 


i(l-J_)V+ 



Darboux I 



ln|q|dq2 J_J+ 



q^ 21nyjr 

_ , l + ln^lql , J-\n^JjZj+ 

Darboux II dq^ ^ ^ 



Darboux Ilia "^4*^' dq' 



q^ln'lql 2{l+\r^^) 



q* 2(1 + ^) 

2(fc + J_) 



Darboux Illb (fc + q^ ) dq^ 



Darboux IV « + cos(ln|q|) , (^'^ 

q2sin2(ln|q|) 2(a + cos(ln i/jT)) 



Taub-NUT 


(4m 


+ |q|)dq2 






|q| 


2 (4m + a/jT) 


!/-fold a ^ 

i/-fold a = 


(a + 


&|q|^)dq2 


i_ j_ 


|q|- 


1 |2— 1 

|qr " 

"2dq2 


2(a + 6J^) 



1)k). The generic coordinates q are identified with Poincare coordinates [41] coming 
from a stereographic projection in R^^-'^. Note that r = |q| is not a geodesic distance; 
in fact, the proper geodesic coordinate f is related to r by 

r = ^ tan ds' = df' + ^j^fi^ d^'. (3.23) 

^/k \ 2 J K 

As it is well-known the geodesic motion on these spaces is MS and their additional 
quadratic integral of motion in Poincare coordinates can be found in [2] . In table 1 
we write the resulting expressions for k — {±1, 0}. 

• ND Darboux spaces. The four 2D Darboux spaces are surfaces of nonconstant curva- 
ture whose geodesic motion is quadratically MS [42, 43]. From the point of view of 
superintegrability, these surfaces can be seen as the closest analogues of the Rieman- 
nian surfaces of constant curvature. An ND spherically symmetric generalization 
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for them was recently proposed in [37] by requiring them to be endowed with an 
sl{2, M)-coalgebra symmetry. We remark that the Darboux surface of type III ad- 
mits two (equivalent) generalizations of this kind: the type Ilia given in [37] and 
the type Illb constructed in [26]. Only for this Darboux III space the MS prop- 
erty has been proven in arbitrary dimension [26] (and the additional integral is also 
quadratic), while for the three remaining types the MS problem is still open. Recall 
also that the parameters a and k appearing in the types Illb and IV in table 1 are 
real constants. 

• Iwai-Katayama spaces. These are the A^D counterpart of the 3D spaces underlying 
the so called "multifold Kepler" systems introduced by these authors in [44]. The 
generic i^-fold Kepler space depends on two real constants, a and 6, as well on a 
ra^ionaZ parameter v as shown in table 1. For the sake of a further more transparent 
discussion we have distinguished the cases a 7^ and a = 0, and in the latter case h 
can always be taken equal to 1. The Iwai-Katayama systems are of physical interest 
as they are generahzations of the MIC-Kepler and the Taub-NUT ones. As far as 
the underlying space is concerned the Taub-NUT metric is recovered from the v-lold 
Kepler one provided that z/ = 1, a = 4m and 6=1 (recall that the MIC-Kepler 
space is just the Euclidean one). We stress that in the 3D case all these Hamiltonians 
have been shown to be MS, but the additional integral of motion is not, in general, 
quadratic in the momenta (see [73] for a detailed discussion on the subject). 

It is also worth stressing that in 3D the classical Riemannian spaces, the Darboux 
space of type Illb and the Iwai-Katayama spaces are just particular cases of the so-called 
"Bcrtrand spaces" . The A^D version of the whole family of 3D Bertrand spaces has been 
recently studied in [50, 73] by making use of the same sl(2, M)-coalgebra symmetry, and 
their explicit expressions for the abovementioned particular cases are already included in 
table 1. For the sake of simplicity, we have omitted the generic form of the Bertrand 
metrics in table 1, but their ND generalization only requires to consider the appropriate 
conformal factor coming from [50, 73] and to replace the metric dO? (and its associated 
2D angular momentum L^) on the 2D sphere by the metric (3.8) on the (A^ — 1)D sphere 
(and associated {N — 1)D "angular momentum" (3.16)). 

4 Superintegrable potentials on curved spaces 

At this point, the Hamiltonian (1.1) can be thought of as the kinetic energy (3.15) plus 
a /i^-term and a function W(|q|), and by considering a symplectic realization with non- 
vanishing 6j's. This system can also be expressed in terms of the spherical phase space 
variables (3.9) and (3.13) by taking into account that the generators {J±, J3} turn out to 
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be 

J_ = J3 = rpr, 



N-l ,2 



. 2 , -2t2 , V- S (4.1) 

In this way we find that the Hamiltonian (1.1) can be rewritten as 

^ = ^^^^^ ^ ^7(0^ ^ ^^^^^ 

^ V + — - U 2) 



Again, we stress that the symplcctic rcahzation of the generator J+ is the essential tool 
to incorporate in a natural way up to N centrifugal potentials associated to the 6j's. 

Therefore this Hamiltonian is QMS as it is endowed with the {2N — 3) integrals of 
motion given by (2.2). However, the existence of an additional integral providing the MS 
property for Ti, is, in general, not guaranteed. Nevertheless, as it already happened with 
the free motion, Ti is separable in the latter form (4.2). In this case the integrals of motion 
C(rn) (2.2), with arbitrary centrifugal terms, are given by 

Consequently, we obtain again a common set of A^" — 1 angular equations 

(a \ 2 , ^N-i , 
C(2)(6^iV-l,Pe^_J = Pg^_^ + + -^-2^ , 

cos t7jv-i sm (7JV-1 

C^i)i9N-i+i,Pe^.,^,) = Pe^_,^, + —2^ + ' I = 3,..., N-l, 

sm ttjsf-i+i cos 

2 Qiv-i) bi 



C(iv)(^i,P.J=pl + -^ + 



sin 61 cos^ 9i 

^ L2 , , /4 4^ 

^cos^e,ni=iSm^^. nl-.'sm'e,' ^'^ 

(to be compared with (3.22)) plus one radial equation which depends on the specific 
Hamiltonian under consideration 

^(.p.)= ^-+'-"g:;r'-""' ^»w. (4.5) 

Thus when any hi ^ the spherical symmetry is broken and the "total angular momen- 
tum" (3.16) is no longer a constant of motion, its role being now played by C(jv) (4.4) 
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(only when all bi = 0, C(^n) reduces to L^). In the case with 6j ^ 0, the generator J+ (4.1) 
simply reads 

J+=p2 + r-2C(jv), (4.6) 
and the Poisson brackets among the s[(2,R) generators (4.1) and turn out to be 

{J_,L2} = {J3,L2} = 0, 

r 7 T 21 _ 4 sin^^tangjpg^- 

X ( 6, tan^^,- - y p. ^ ^ I . (4.7) 

,4^,cos2^,ni-4+i«in^^. nf=-+isin^^J 

Summing up, once a single nonzero parameter 6j is allowed, the spherical symmetry is 
broken, and the chain of rotation subalgebras so{m) C so{N) do not provide symmetries 
for the Hamiltonian, as it was the case of the Casimirs (3.19) for the free motion. There- 
fore, the s[(2, M) Lie-Poisson coalgebra symmetry can be understood as the appropriate 
generalization of the usual rotational symmetry which allows to deal with the general case 
above. 

We also remark, as pointed out in the introduction, that the yU^-potential in (4.5) (or 
(4.2)) can be interpreted as a Dirac monopolc-typc term [44, 56, 57, 60]. Hence since 
C(Ar) is an integral of the motion, the appearance of this monopole-type potential can 
be seen, from a dynamical viewpoint, as coming from a radial centrifugal term ruled 
by C(iv)- In particular, when there are no centrifugal potentials {b^ — 0), the /x^-term 
has already been interpreted, for some concrete systems, as a proper Dirac monopole 
[44, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61] that actually comes from the "total angular 
momentum" through C(Ar) = — >• + /i^. 



5 KC and oscillator potentials on curved spaces 

Among all the possible choices of the central function W(r), two of them should correspond 
to the appropriate iVD definition of the intrinsic KC and oscillator potentials on Ai. As 
we shall see in the sequel, a solution to this problem can be found through a suitable 
generalization of the 3D construction for such potentials. 

Let us consider the Laplace-Beltrami operator on the 3D spherically symmetric man- 
ifold A4^ with metric (3.7) in the coordinates q: 

where g''^ is the inverse of the metric tensor Qij and g is the determinant of Qij. The 
radial symmetric Green function t/(|q|) = U{r) on JA^ (up to multiplicative and additive 
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constants) is defined as the positive nonconstant solution to the equation 



A^,3^(r) = ^A(^rV(r)^)=0 on M^O} , (5.9) 

so that 

^(^)-f\^y ("») 

Now, the essential point in the A^D case [26] is the fact that we can keep the very same 
definitions given in [46, 47, 48, 49, 50] for the KC and oscillator potentials on a 3D 
spherically symmetric space. In particular, the intrinsic KC potential on the A'^D space 
A4 will be defined by 

WKc(r) -.^aUir), (5.11) 

while the intrinsic oscillator potential is defined to be proportional to the inverse square 
of the KC potential 

^o(r) := jj^y (5.12) 

where a and f3 are real constants. 

Clearly, the intrinsic KC and oscillator potentials do preserve the s[(2, M)-coalgebra 
symmetry since they are defined through a function of |q| = r = y/JZ. As a straight- 
forward consequence, such a symmetry allows us to give the superposition of either the 
intrinsic KC (5.11) or oscillator (5.12) potential with A^ centrifugal terms and also with 
a monopole-type term. In terms of the symplectic realization (2.3) and by considering 
arbitrary parameters bi, the curved KC system read 

2/(v^)2 +^Kc(v^-) - — 27cra7 — + "^^1^1^ ^^-^3) 

and the intrinsic oscillator on a curved space is defined as 

2f(vrjf +"°(^-'-' vMf *uMf- 

Both systems can be immediately written in terms of the spherical phase space variables 
by means of (4.1). 



6 Some physically relevant examples 

The results of the previous section are exphcitly illustrated by writing in table 2 the QMS 
Hamiltonians (5.13) and (5.14) associated to all the particular spaces given in table 1. 
For each A^D curved space we display the KC, oscillator, monopole-type and centrifugal 
potentials (with parameters a, (3, fi^ and 6,, respectively). A specific QMS Hamiltonian 
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Ti thus comes from the superposition of some (or all) of the corresponding terms. Notice 
that, in some cases, the multiplicative constants that arise when computing the Green 
function (5.10) have been absorbed within a and (5. 

In what follows we explicitly analyse these systems for the three families of spaces 
discussed in section 3.2, since they are indeed the most relevant ones from the physical 
point of view. 



Table 2: Intrinsic KG, oscillator, monopole-type and centrifugal potentials with sl(2, M)- 
coalgebra symmetry in the generic coordinates q corresponding to the free Hamiltonians 
on the particular A'"D spherically symmetric spaces given in table 1. 



Space 


KC 


Oscillator 


Monopole 


Centrifugal terms 


Euclidean 

Spherical 

Hyperbolic 


a 

~H 

a(q2 - 1) 




2q2 

M'(l + q2)2 


2—1 ^ 

i=i ^» 
1 ^ h 

i(i-<.')T| 

1=1 * 


|q| 

a(q' + 1) 


(q2 - 1)2 
^q' 


2q2 




(q2 + 1)2 


2q2 



(3 ^ q^ ^ 6 

Darboux I a\/hi Iql - — — , , , 



ln|q| 21n|q| 21n|q|^g2 

n K TT /i , , 2, I /3 M^in^|q| q^in^|q| f &i 

Darboux 11 a\/ 1 + m q ^ o o > -o 

V l + ln'|q| 2(l + ln2|q|) 2(1 + In^ |q|) £^ g,^ 



Darboux Ilia a\/l + |q| 



P M^q^ q^ f h_ 

l + |q| 2(l + |q|) 2(l + |q|)^g2 



Darboux Illb "^"^ 



4- n2\Z^ „2 



|q| A; + q2 2q2(fc + q2) 2(fc + q2)^g2 

u n I /3 /x2sin2(ln|q|) q2 sin2(ln |q|) f 6i 

Darboux IV aV <i + cos(ln q ) — --- — . inn ttt , ,,, > 

^ ^ ' a + cos(ln|q|) 2(a + cos(ln |q|)) 2(a + cos(ln |q|)) ^ g2 

Taub-NUT av/4m|q|-i + l „ , „ :^7^^^^^ V ^ 

^ ' ' 4m+|q| 2|q|(4m+|q|) 2(4m + |q|) ^ ^2 

i^-iold a a\ a\a\ " + o t— > 

^ a|q|-^+6 2|q|^(a + 6|q|i) 2(o + 6|q|i) ^ 9? 

Iql" 2|q|. 2 
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6.1 The classical Riemannian spaces of constant curvature 



For any simply connected Riemannian space with constant sectional curvature the 
metric function / in Poincare coordinates [2] reads 

/(|q|) 



1 + Kq2' 

and the Green function (5.10) is easily computed: 

/tq2 — 1 



t^(|q| 



|q| 



As a consequence, the KG and oscillator potentials are the ones shown in tabic 2 for 
K = {±1, 0}. Such expressions can be rewritten in a more usual form by introducing the 
geodesic radial coordinate r (3.23) which gives 



7j ^ 11 ('■\ fltan2(^f) 



tan(y^f) ' K 

We recall that when all the iV parameters hi are diferent from zero, the generalized KG 
Hamiltonian is QMS but not quadratically MS. Moreover, only when at least one of the 
6j's vanishes an additional quadratic integral of motion arises [2] (which is a component of 
the Laplace-Runge-Lenz A?"- vector). Nevertheless, it has been recently proven in [36] that 
when the KG potential is constructed on the 3D Euclidean space it is possible to consider 
the three centrifugal potentials obtaining a MS Hamiltonian with an additional quartic 
integral. In contrast, the so-called Smorodinsky-Winternitz system [i.e., the oscillator 
plus centrifugal terms) is always MS for any value of the 6i's. 

We stress that the MIG-Kepler Hamiltonian [44, 56, 60, 62, 63, 64] is also recovered 
when the KG and the monopole potential are considered in the fiat Euclidean space: 

Euchdean : T^Mic-Kepicr = -J+ 7^ + 7rT~ = oP ^ Tl + (^-1^) 

z Y iV— zj_ z |q| zc^ 



Therefore, as a byproduct of our procedure, the curved MIG-Kepler analogue on the A"D 
spherical and hyperbolic spaces can be constructed from table 2. Namely, 

Spherical : ^Mic-Kepier = ^(1 + J-fJ+ + '^^^^ ' ^ ^'"^ 



2' ' ^ VJ_ 2J_ 

ov |q| 2q^ 



1 V n, 1/. T ^2 X a(^- + l) ^x\l-J-f 

Hyperbohc : TtMic-Kepier = -(1 - J-) J+ -j== \ — 
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It can be easily checked that in this way we have exactly recovered (in Poincare coor- 
dinates) the curved MIC-Kepler systems studied in [74, 75, 76]. Obviously, all these 
Hamiltonians can be generalized by adding the 6i-centrifugal terms given in table 2, and 
in that case the QMS property is, by construction, fully preserved. 



6.2 Darboux spaces 

A thorough discussion of the MS potentials on the four types (1, 11, llla and IV) of 2D 
Darboux surfaces was given in [43]. Hence it is natural to compare with them the 2D 
versions of the A'"D potentials given in table 2. To carry out this analysis, we firstly 
recall that the Darboux metrics given in [43] depend on two coordinates {u, v), and their 
associated metrics and free Hamiltonians read 

ds^^F{uf{du^ + dv^), H^^;^, (6.18) 

where {pu,Pv) are the conjugate momenta and the function Fiu)"^ is given by 
Type I: F{uf = u; Type II: F{uf = 1 + -u"^; 



Type Ilia: ^(m)^ = e"^" + e""; Type IV: ^(m)^ = 

sin u 



(6.19) 



Secondly, we also recall that the A^D spaces given in table 2 are just an NT) spherically 
symmetric generalization of these four spaces that was constructed in [37] through the 
maps 

M^lnr = ln|q|, dv'^^d^f. (6.20) 

Now, if we consider our A^D KC and oscillator potentials given in tabic 2 for these 
four spaces and we perform the substitution -u = In |q|, we immediately sec that, in any 
dimension, we recover expressions for the potentials that depend only on the variable u. 
Now, if we go back to the classification given in [43], the unique MS potentials in 2D 
which are functions of u alone are 



Type I: -^VTl\ Type II: ^ — -L^; 

u ln|q| 1 + 1 + ln |q| 

11 11 

Type Ilia: > j — r; Type IV: 



(6.21) 



l + e" l + lql' ' a + cosii a + cos(ln |q| 



Surprisingly, we find that these four potentials are just the intrinsic oscillators of 
the four 2D Darboux spaces, and there is no KC potential from table 2 appearing in the 
classification given in [43]. This fact suggests that in the case of 2D spaces of non-constant 
curvature, the intrinsic oscillator potential seems to be more fundamental than the KC 
one from the integrability viewpoint, since the former would be a MS system whilst the 
later would be only a QMS one. 
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A remark is in order: in the classification [43], the MS potentials associated to the 
Darboux spaces II and Ilia contain also other terms depending on u, namely 1/(1 + u~'^) 
and 1/(1 + e~"), respectively, which in fact can be obtained by adding a constant 7 to 
the above oscillator potentials with constant /3: 

P , /5 + 7 , 7 



T TTT , /^ + 7 , 7 ^ ^ 

lype 111a: h 7 = 1 . 

1 + e" 1 + e" 1 + e-" 

Finally, recall that in arbitrary dimension A^, only the oscillator potential with N cen- 
trifugal terms for the Darboux space lllb has been proven to be MS [26]. Therefore the 
MS property remains as an open problem for all the remaining types of ND intrinsic 
oscillator potentials. 



6.3 Iwai— Katayama spaces 

We have also written in table 2 the resulting potentials corresponding to the Iwai- 
Katayama spaces described in section 3.2 by considering either a 7^ or a = (with 
6=1). We remark that, although the z^-fold Kepler metric given in table 1 depends con- 
tinuously on a, a glance at the potentials (displayed in table 2) reveals that it is convenient 
to perform a separate analysis of the cases a and a — 0. 

Firstly, let us recall that the A^D version of the 3D multifold Kepler Hamiltonian 
introduced in [44] is given by 

7^.foid-Kepier = ^ fp' + A*^q"' + A*'c \q\^-^ + f/d |q|^-') , (6.23) 

2(a -I- olql") ^ ^ 

where z/ is a rational parameter and a, b, c and d are real constants. By expanding this 
expression we find that 

/t,yfold-Kepler — — ,i, + ^, , , i — + ^, , i , .1, + — , . . 1 , ■ (,0.24j 

2{a + b\q\T:) 2{a\q\-T:+b) 2|q|i;(a + 6|q|i;) 2(a + 6|q|^^) 

Hence, table 2 allows us to provide a more clear interpretation of the four terms forming 
this Hamiltonian, which indeed does depend on the value of the constant a. Namely: 

• If a 7^ 0, the first term in the Hamiltonian (6.24) is the kinetic term written in 
table 1, the second is an intrinsic oscillator with (3 = fi'^d/2, the third is the Dirac 
monopole and the fourth comes out by adding a constant 7 to the corresponding 
oscillator potential (this trick is the same as the one used in (6.22)): 

/5 . /5 + 67 07 

+ 7 = ^ 7 + TT^^ (6-25) 



a|q| '^+b a\q\ + b a-|-6|q| 
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so that 7 = n'^c/{2a). Consequently, from this approach we can state that the 
multifold Kepler systems with a ^ are, in fact, multifold oscillator systems. This 
interpretation was already given in [50] for the 3D case. 

• If a = (and 6=1), the multifold Kepler Hamiltonian (6.24) reduces to 

1 2 A*^*^ A*^ A*^^ 

ni/fold-Kepler; a=0 = 7^|q| P + -— T + — + -pT ■ [6.26) 

^ 2\q\i' 2\q\'^ ^ 

Therefore in this case the first term is the kinetic energy given in table 1, the second 
is an intrinsic KC potential with a — — /x^c/2, the third is the monopole and the 
fourth is an additive constant. Thus only in the case a — the Hamiltonian (6.24) 
is a proper multifold Kepler system. 



From the viewpoint adopted in this paper, it is apparent that the case a 7^ corre- 
sponds to the u-iold generahzation of the Taub-NUT system [44, 59, 65, 66, 67, 68, 69, 70, 
71, 72], while a = is the ly-Md version of the MIC-Kepler model [44, 56, 60, 62, 63, 64]. 
Explicitly: 



The Taub-NUT system arises as the particular case of (6.24) with v — 1, a — 4m, 
b = 1, c= l/(2m) and d = l/(4m)^ (see [44]) which yields 

^ _ |q|p^ , /i^|q|/(4m)^ ^2 ^ ;i7(4m) 



2(4m + |q|) 2(4m+|q|) 2|q|(4m+|q|) 4m+|q| 

+ (6.27) 



2(l + 4m/|q|) 2(4m)2 V |q 

so that the same interpretation as for the multifold Kepler system with a ^ holds. 
In the first line of (6.27), the first term is the kinetic energy, the second is an 
intrinsic oscillator with f3 = yU^/(2(4m)^), the third is the monopole and the fourth 
corresponds to adding a constant 7 = /i'^/(4m)^ to the oscillator potential, leading 
to /3 ^ /3 + 7. 

The (fiat) MIC-Kepler system is recovered from (6.24) or (6.26) when 1/ = 1, a = 
d — 0, b — 1 and c = — 2q;//x^, giving rise to the expression (6.15). 



Note also that in 3D, the generic Hamiltonian (6.23) (for any value of a) has been re- 
cently shown to be MS [73] . This Hamiltonian has an additional integral of motion (com- 
ing from a generalized Laplace-Runge-Lenz vector) which, in general, is not quadratic in 
the momenta. 

Let us end by stressing that all the A^D Hamiltonians that wc have presented are 
QMS, but some of them could be MS, as we have already commented. As a matter 
of fact, one of the models that we have described (namely, the Darboux HI oscillator 
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potential with monopole and centrifugal terms) is, to the best of our knowledge, the only 
known quadratically MS system living in an A^D space of nonconstant curvature; the 
additional integral can be consulted in [26]. The search for an additional integral for 
the rest of the systems in the Darboux spaces presented here remains a challenging open 
problem. In the cases where such a constant of the motion exists, it cannot be derived 
from the sl(2, M)-coalgebra symmetry, and it reflects the exceptional superintegrability 
properties of the Riemannian manifold determined by f{r). 

On the other hand, it is worth remarking that the results presented throughout the pa- 
per can be extended as well to curved Lorentzian spaces through an analytic continuation 
method analogous to the one applied in [37] for some free Hamiltonians. 

Finally, we would like to emphasize that the underlying 5((2, R)-coalgebra symmetry 

can also be implemented in the Quantum Mechanical analogues of the systems here pre- 
sented. In particular, the quantum counterparts of the integrals of the motion (2.2) can be 
readily obtained after dealing with the ordering problems which arise in the quantization 
of (1.1) due to the term p^//(|q|)^. We shall report on these and other related issues 
elsewhere. 
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